Spin-lattice model of Magneto-electric Transitions in RbCoBr 3 
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Extensive Monte Carlo simulations are performed to analyze a recent neutron diffraction experi- 
ment on a distorted triangular lattice compound RbCoBr3. We consider a spin-lattice model, where 
both spin and lattice are Ising variables. This model explains well successive magnetic and dielectric 
transitions observed in the experiment. The exchange interaction parameters and the spin-lattice 
coupling are estimated. It is found that the spin-lattice coupling is important to explain the slow 
growth of a ferrimagnetic order. The present simulations were made possible by developing a new 
Monte Carlo algorithm, which accelerates slow Monte Carlo dynamics of quasi-one-dimensional 
frustrated systems. 

PACS numbers: 75.80.+q, 75.40.Mg. 77.80.-e 



I. INTRODUCTION 

Frustrated magnets have been attracting much interest 
for many decades X The ordinary magnetic order is de- 
stroyed and the ground state may remain disordered or 
turn into an exotic state. The ground state of a frustrated 
system is usually unstable against a small perturbation. 
The system manages to find a way to relax frustration 
and change the state. We may design and control mate- 
rial functions using the frustration effects. 

The ABX3-type compounds are well-known frustrated 
magnets. The lattice structure is the stacked triangular 
lattice. There is frustration when the nearest-neighbor 
interactions are antiferromagnetic. Successive magnetic 
phase transitions occur because of the strong frustra- 
tion when the spins have the Ising anisotropyi^ The 
low-temperature magnetic structure is the ferrimagnetic 
state. There exists a partially-disordered (PD) phase 
between the paramagnetic phase and the ferrimagnetic 
phase. In the PD phase, one of the three sublattices 
is completely disordered, while the other two sublattices 
take antiferromagnetic configurations. There is no struc- 
tural phase transition in most compounds. The system 
remains fully frustrated down to the lowest temperature. 

The KNiCl3-family compounds are exceptional in that 
they exhibit structural phase transitionsj 4 ' 5 ' 6 We can 
observe the structural phase transitions by the dielec- 
tric measurements because each BX3 chain has a neg- 
ative charge. These compounds have both magnetic 
and dielectric characteristics: we call them the magneto- 
dielectric compounds. The magnetic phase transitions 
and the structural (dielectric) phase transitions usually 
occur at different temperatures. However, Morishita et 
al. found that both transitions occur at the same tem- 
perature in RbCoBr 3 j£i£ It is a very rare case among 
the KNiCl3-family compounds. Magnetic and dielec- 
tric measurements^^ found that the phase transitions in 
RbCoBr 3 are quite unusual in the following points com- 
pared with other compounds: 

(i) The dielectric transition temperature 37 K of 



RbCoBr3 is very low compared with other com- 
pounds, for which the transition takes place around 
the room temperatures. The energy scale of the 
structural (dielectric) system in RbCoBr3 seems to 
be suppressed somehow. 

(ii) The temperature dependence of the dielectric con- 
stant does not exhibit a diverging behavior. This is 
clearly different from another KNiCl3-family com- 
pound RbFcBr3, which exhibits sharp divergence 
at the transition temperature, 34.4 K£ 

(iii) The increase of the spontaneous polarization below 
the dielectric transition temperature is very slow, 
while that of RbFeBr3 is very sharp. 

(iv) The magnetic PD phase appears in a very narrow 
temperature region. The first neutron measure- 
ment suggested that it might disappear^ A recent 
improved neutron experiment^ 1 - made it clear that 
it exists between 31 K and 37 K. This is also a clear 
difference from other compounds such as CsCoBr^ 
and CsCoCl 3 4 

(v) The growth of the ferrimagnetic order is very slow. 
The neutron count increases linearly with the tem- 
perature decrease in the low-temperature phase. 

These characteristic behaviors suggest that there is an 
unknown mechanism of interplay between the magnetic 
system and the dielectric system in RbCoBr 3 . 

An aim of this paper is to propose a proper theoret- 
ical model that quantitatively explains the experimental 
results of RbCoBr3. A well-known theoretical model for 
the ABX3 compounds is the antiferromagnetic spin sys- 
tem on the stacked triangular lattice. 12 : 13 : 14 ! 15 ! 16 How- 
ever, the ordinary spin model without a coupling to the 
dielectric system is not sufficient to explain RbCoBi'3. 
The chain-mean-field theory^ gives the magnitude of the 
second-nearest-neighbor magnetic interactions on the ex- 
plain IJ2I ~ 1 K, which is comparable to that of the 
nearest-neighbor interactions |Ji| — 2.5 K£ This is not 
acceptable from the experimental point of view. 
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In the present paper, we use the spin-lattice model pro- 
posed by Shirahata and Nakamura>il This model showed 
that the PD phase may disappear because of the relax- 
ation of frustration by the lattice distortion; each of the 
spin system and the lattice system relaxes frustration of 
the other. A single transition may occur from a param- 
agnetic and paraelectric phase to the ground state phase 
without experiencing the intermediate PD phase. Shi- 
rahata and Nakamura also noticed that the cooperation 
between the spin system and the lattice system works 
only when the energy scale of the lattice system is com- 
parable with that of the spin system. 

In this paper, we refine the above spin-lattice model. 
We find that a soft lattice system coupled with a spin 
system explains the interesting material RbCoBr 3 . We 
thereby clarify the origin of the characteristic behaviors 
of this compound. For the purpose, we develop a new 
Monte Carlo (MC) algorithm, which eliminates slow MC 
dynamics in quasi-one-dimensional frustrated spin sys- 
tems. We have performed extensive MC simulations and 
determined various physical parameters. 

We explain our model Hamiltonian in Sec. [TTJ A nu- 
merical method is explained in Sec. IIII1 and the results 
are presented in Sec. IIVI Discussions are given in Sec. 

El 



II. THEORETICAL MODEL 
A. Structure of ABX3 compounds 

The lattice structure of ABX3-type compounds is the 
stacked triangular lattice. Face-sharing BXg octahedra 
run along the c-axis forming a BX3 chain. Magnetic 
B 2+ ions form an equilateral triangular lattice on the 
c-plane, which causes frustration. Exchange interactions 
along the BX3 chains, J c , are antiferromagnetic. The 
magnitude of J c is much larger than that of the nearest- 
neighbor interactions J\ on the c-planc: \J C \ \Ji\. 
Therefore, this spin system can be considered as a quasi- 
one-dimcnsional system with frustration on the c-planc. 

The typical lattice structure of ABX3-type compounds 
at high temperatures is shown in Fig. [1] (a). The space 
group is POt^/mmc. Magnetic ions forming an equilat- 
eral triangular lattice sit on a level plane. This struc- 
ture remains down to the lowest temperature in most 
compounds. In the KN1CI3 family, the structural phase 
transitions occur as we decrease the temperature. Each 
BX3 chain shifts upward or downward keeping the rel- 
ative distance between the atoms. One of the lattice 
structures after the structural phase transitions is shown 
in Fig. [T] (b), where two sublattices on the triangular 
lattice shift upward with the same amount while one 
sublattice shifts downward. The space group is P63CTO. 
It is the ferrielectric structure of KMCI3 observed at 
the room temperatures^ We refer to this structure as 
"lattice-Ferri" or "f-T-J." m this paper. Another possible 
structure is a configuration with one sublattice shifting 




FIG. 1: Typical crystal structures of ABX3-type compounds. 
A-ions are omitted. Black circles depict magnetic B 2+ ions 
and grey circles depict X~ ions. Each arrow depicts the 
shift direction of a chain, (a) A symmetric structure at high 
temperatures. The space group is P63 /mine, (b) A room- 
temperature KMCI3 structure. The space group is P&zcm. 
We call this structure "lattice-Ferri". (c) A low-temperature 
structure. The space group is P3cl. We call this structure 
"lattice-PD" . (d) Another low-temperature structure. The 
space group is P3cl. We call this structure a "three-sublattice 
lattice-Ferri". Each sublattice polarization takes a different 
value. 

upward, one shifting downward and the third unchanged 
as shown in Fig. Q] (c). The space group is P3cl. We 
refer to this structure as "lattice-PD" or "t-|-0" in this 
paper. Nishiwaki and Todoroki^ discussed the appear- 
ance of the three-sublattice ferrielectric state in RbCoBr3 
using the mean-field approximation. It is a structure with 
the t"T"l configuration but the amount of displacement 
in each sublattice is different from the others, as shown 
in Fig. Q] (d). The space group is P3cl. We refer to this 
structure as "three-sublattice lattice-Ferri" in order to 
distinguish from the "lattice-Ferri" structure ofFig.fjJb), 
which has the two-sublattice order. 



B. Spin-lattice model 

We consider a model on the stacked triangular lattice 
with spin and lattice degrees of freedom^ The size in 
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the a and b directions is J, while the size in the c direc- 
tion is L c . There are a spin variable iSy and a lattice 
variable at each site. Here, the subscript i denotes 
the position in the c-axis, while j denotes the position 
on the c-plane. We define each spin as an Ising variable 
Sij = ±1/2 because Co ions have the Ising anisotropy. 

The lattice variable denotes the displacement from 
the symmetric lattice point along the c-axis. We ap- 
proximate the displacement with the Ising variable as 
(Tij = ±1/2; each ion shifts either upward (er^ = 1/2) 
or downward (cr.y = —1/2). The reason of approximat- 
ing the present lattice system with the Ising variables 
is as follows. The symmetric lattice structure with zero 
displacement as shown in Fig. QJa) (Pds/mmc) appears 
at high temperatures. Structural phase transitions occur 
successively as the temperature decreases. The lattice 
configuration is the [-|-0 state (Fig. [He)) m the inter- 
mediate phase and is the [-[-[ state (Fig. [TJb)) in the 
low-temperature phase. This is analogous to the succes- 
sive magnetic phase transitions from the paramagnetic 
phase to the PD phase and to the ferrimagnctic phase in 
the Ising model on the stacked triangular lattice. In the 
previous paper p£ we considered a lattice variable taking 
three states, +1,0,-1. Here, we omit a state (Xjj = 0. 
The (Jij = state can be represented by a mixture of the 
Uij = 1/2 state and the ov, = — 1/2 state. A chain shift 
is the sum of the ion displacements along the chain in the 
experiment. 

The Hamiltonian consists of the lattice part 7Yl and 
the spin part 7is : 

H = H L +H S , (1) 

where 

n.n. 

W L = - 2J c L ^cr i:) cr( i+1)j ■- 2Ji^2^2<7 ij a lk 
hi i <jk) 

n.n.n. 

1 (jk) 

Hs = - 2jS^2SijS {i+1)j 

hi 

n.n. 

- 2jf £(1 - A(ofc - a lk f)S l3 S lk 

n.n.n. 

- 2J IJ2 Yl t 1 - - VikftSijSik- (3) 

The lattice part comes from the elastic energy: (cr.y — 
ai'j') 2 . The spring constant is denoted by l 2 y where 
each of (c, 1, 2) denotes a direction of the interaction: 
c denoting along the c-axis, 1 denoting the nearest- 
neighbor(n.n) pairs on the c-plane, and 2 denoting the 
next-nearest-neighbor(n.n.n) one on the c-plane. The 
sign of the spring constant is determined based on the 
effect of the exclusion volume effect. It is positive along 
the c-axis: J r L > 0. An ion pushes the next ion in the 



same direction. The spring constant for the nearest pairs 
in the c-plane should be negative: j\ < 0. An ion shifts 
upward if the neighboring ion shifts downward, because 
ions try to stay away from the neighboring ions. There- 
fore, there is frustration in the triangular lattice. Wc 
choose J 2 to be positive in order to realize the [-[-[ 
state observed experimentally at low temperatures. 

The main idea of this paper is to make the spin-spin 
exchange integrals depend on the lattice variables. We 
assume that the interaction becomes weak if the exchange 
path is distorted. Thus, the in-plane exchange interac- 
tion becomes (1 — A) times weaker if aij and <ii k have 
opposite signs. We use the same value of A for jf and 
jf for simplicity. The exchange path along the c-axis is 
rigid against the ion shift and hence we assume Jf to be 
unaffected. 

We assume that the nearest-neighbor spin-spin interac- 
tion is antiferromagnetic ( Jf < 0) and the next-nearest- 
neighbor interaction is ferromagnetic ( Jf > 0) in order to 
realize the ferrimagnctic state in the ground state. The 
interactions along the c-axis in the real compound are 
antiferromagnetic. (jf < 0) 

The lattice part and the spin part of the Hamiltonian 
have the same form of the antiferromagnetic Ising model 
on the stacked triangular lattice. They are connected 
by the A term of the form — 4jf 2 ^ <J ij' J ikSijSik. Thus, 
the present model can be regarded as the Ashkin- Teller 
modelJ^ 

Experimental estimates of the exchange integrals were 
J c s ~ -62 K, Jf ~ -2.5 K, and jf ~ 1 K4 The estimate 
of Jf was obtained from the position of the broad maxi- 
mum peak of x\\- We consider it underestimated, which 
we will discuss in Sec. HVl 



C. Relaxation of frustration by lattice distortion 

Here, we consider how the ordered magnetic state is fa- 
vored by the lattice distortion in the present spin-lattice 
model. When the lattice takes the lattice-Ferri (T-?-|) 
configuration, the PD state of the spin system is fa- 
vored magnetically As shown in Fig. [Ha), the nearest- 
neighbor interactions between two f-shifted sublatticcs 
remain strong (depicted in the figure by thick lines), while 
those between an | -shifted sublattice and a |-shifted sub- 
lattice are weakened by the A term. The strong bonds 
form a honeycomb lattice, where the spins are ordered an- 
tiferromagnetically. The remaining spins on the .[-shifted 
sublattice interact with the spins on the honeycomb lat- 
tice through weak bonds. The molecular field on the 
[-shifted sublattice from the [-shifted sublatticcs van- 
ishes because of the antiferromagnetic ordering on the 
[-shifted sublatticcs. Then, the spins on the [-shifted 
sublattice may be disordered. We call the PD state of 
the spin system as "spin-PD" in this paper. 

The similar argument is possible when the lattice 
system takes the lattice-PD ([-|-0) state as shown in 
Fig- Ob). The ferrimagnctic spin state is favored in this 
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FIG. 2: (Color online) Relaxation of frustration by the lattice 
distortion. Red (solid) circles depict c-chains shifting upward. 
Blue (grey) circles depict c-chains shifting downward. Open 
(white) circles depict c-chains not shifting. Thin (thick) lines 
depict weak (strong) interactions, (a) When the lattice is 
deformed to a lattice-Ferri (j-j-J.) pattern, the spin-PD state 
is favored, (b) When the lattice is deformed to a lattice-PD 
(T-J.-0) pattern, the spin-Ferri state is favored. 

case. The nearest-neighbor interactions between an |(4) _ 
shifted sublattice and a 0-shifted sublattice are weakened 
by A/4, while those between the j-shifted sublattice and 
the ^-shifted sublattice are weakened by A. The anti- 
ferromagnetic ordering is realized on the stronger bonds, 
which is the ferrimagnetic state. We call it in this paper 
"spin-Ferri" . The present mechanism was discussed by 
Plumer et al.r^ when the lattice distortion is static. 



III. MONTE CARLO METHOD 
A. Axial-bond-cluster flip algorithm 

We briefly explain our new Monte Carlo simulation 
algorithm. The detail will be reported elsewhere^ 

The origin of the slow MC dynamics in the quasi-one- 
dimensional ( | J c ^> | Ji | ) Ising system is a very long 
correlation length along the c-axis. It rapidly grows at 
low temperatures as £ c ~ exp[| J c \/T}. A large magnetic 
domain along the c-axis is not flipped by the standard 
single-spin-flip algorithm. Koseki and Matsubara^ in- 
troduced a cluster-heat-bath algorithm in order to solve 
this problem, but it costs a long CPU time. The pos- 
sible size of simulations is restricted to \J C /Ji\ = 10, 
N = 36 x 36 x 360, and 2 x 10 6 MC steps^ 

Here, we solve this problem using the loop algorithm 
of quantum Monte Carlo (QMC) simulations^ In the 
QMC simulation, a d-dimensional quantum spin system 
is mapped to a (d+l)-dimcnsional classical spin system^ 



i i i 



FIG. 3: A schematic diagram of the MC updating procedure. 
A cluster in an i-chain is defined between two edge arrows. It 
is flipped using the sum of molecular fields from other chains. 



before actual simulations. The additional dimension is 
called the Trotter direction. Then, the classical spin 
system for the QMC algorithm can be interpreted as a 
stacked Ising model. The Trotter direction of the QMC 
system is now the c-axis of the stacked Ising model, and 
the real-space directions of the QMC system are the c- 
plane of the stacked Ising model. The loop algorithm 
of the QMC simulation^ is to flip a "loop" , or an ax- 
ial aligned-spin cluster along the Trotter direction in the 
QMC system. Therefore, the algorithm can be readily 
applied to flip a correlated spin cluster along the c-axis 
of the stacked Ising system. 

The size of the cluster is a stochastic variable in each 
update of the cluster algorithm. Using a proper prob- 
ability we generate locations of the cluster edges, which 
are memorized in the computer array. Then, we calculate 
the sum of molecular fields from other spins to the clus- 
ter between two neighboring edges (two arrows in Fig. [3]). 
The cluster is flipped using the heat-bath probability by 
this molecular field. 

We noticed that the required computational proce- 
dures and the amount of computer memory are inde- 
pendent of the correlation length £ c . We do not need to 
memorize spin states of all sites. Locations of the cluster 
edges and the spin state at each chain edge are stored 
and utilized in the simulation. The linear size along the 
c-axis, L c , is set to £ c times larger than L. The sys- 
tem with L x £ C _L spins is simulated with an effort of 
L 3 . The new algorithm becomes exponentially efficient 
at low temperatures. In this paper, we set L = 104 for 
all data. An effective spin number at low temperatures 
exceeds 10 8 . The periodic boundary conditions are im- 
posed on the lattice. 



5 



B. Observables 

We observe in the present MC simulations the follow- 
ing physical quantities: the sublattice order parameters, 
1/3-structure factors, 1-structurc factors, and the uni- 
form magnetic susceptibility. The sublattice order pa- 
rameters are the sublattice polarization, mh, and the sub- 
lattice magnetization, m 3 , respectively. They are defined 
as 

iV sub 
iv sub . 

where 77 = a, (3, 7 denotes one of three sublattices in the 
triangular lattice, and iV su b = N/3. 

The following structure factors are defined in order to 
detect phase transitions and to compare with the neutron 
experimental data: 




(fi) 2 = <K + m^ + r^) 2 >, (8) 

(A 3 ) 2 = ({m*+m$ + n$)*). (9) 

The 1/3-structure factor takes a finite value when the 
ferrimagnetic state or the PD state is realized. It detects 
the phase transition between the PD phase and the para- 
magnetic phase. The phase transition between the PD 
phase and the ferrimagnetic phase is detected by f\. 

C. Mean- field- like treatment of MC update 

In the present simulation, spin variables and lattice 
variables arc updated separately and alternatively. In 
the calculation of the heat-bath update probability, we 
use the following approximation to simplify the simu- 
lation. For an update of a spin variable Sij, we cal- 
culated the four-body energy, —Ajf 2 &- (J ij (J ikSijSik by 
replacing the lattice variable cry with a mean value 

<fj = Yli=i a ij/Lc, namely as -4Jf 2 A(<fj(fk)SijSik. For 
an update of a lattice variable <7y, we replaced the spin 

variable with a mean value Sj = Sij/L c , namely as 

This mean-filed treatment may be justified by the fol- 
lowing argument. In a cluster updating procedure, the 
sum of the molecular field to a cluster is calculated to 
estimate the updating probability. Since the cluster size, 
~ exp[| jS'Lj/T], is very large around/below the critical 
temperature, the mean over a cluster can be approxi- 
mated by the mean over the whole chain. 



The above mean-field treatment possibly influences the 
critical properties of the phase transitions. Since our 
main purpose here is to explain the experimental re- 
sults, most of the simulations are carried out in the off- 
critical regions. Therefore, we consider that this mean- 
field treatment does not affect our numerical results in 
the present paper. The investigations on the critical 
properties are left for future study. 



D. Simulation conditions 

The choice of the initial state is important in the 
present simulation. Since each of the lattice system and 
the spin system exhibit two successive transitions, we 
have several possible combinations of ordering patterns 
as we change the temperature. We therefore used the 
mixed phase initialization ; 26 ' 27 ! 28 where we prepare sev- 
eral initial spin-lattice states and spatially mix them. For 
example, we start the simulation with the following initial 
state when the temperature is near the spin-PD transi- 
tion temperature. The lattice system above this tempera- 
ture takes the lattice-PD (T-J.-0) state. Because the spin- 
PD state favors the lattice-Ferri (T-t-|) state (Fig.OJa)), 
they may appear at the same temperature. Therefore, 
a half of the system is set to the spin-PD state and the 
lattice-Ferri state, while the other half is set to the spin- 
paramagnetic state and the latticc-PD state. The former 
one appears below the transition temperature, while the 
latter appears above it. We tried other choices of mixed 
states and verified the equilibration. 

The typical number of initialization MC steps was one 
thousand and that of total MC steps was ten thousands. 
It is sufficient except for the vicinity of the transition 
temperature. We performed thirty independent MC runs 
and took the average over these runs. 



IV. RESULTS 

A. Requirements from the experiments 

Experimental findings are listed in the following. They 
should be reproduced by the simulations. 

(i) The uniform magnetic susceptibility shows a broad 
peak at T = 100 KJ£ 

(ii) The dielectric constant shows a small anomaly at 
T = 90 K,Li2 where the latticc-PD (T-j-0) state is 
considered to appear. 

(iii) As the temperature decreases from the room tem- 
perature, the first magnetic phase transition occurs 
at Tni = 37 K. The neutron-scattering data of (1/3 
1/3 1), which corresponds to (/f/ 3 ) 2 of Eq. (J7J), 
show a rapid increase below this temperature, while 
those of (1 1 1), which corresponds to (/f) 2 of 
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FIG. 4: A phase diagram of RbCoBr3. Solid lines show the 
known phase transition temperatures. A broken line at 32 K 
is a new structural transition found in this paper. 



Eq. ©, remains zero. The spin-PD state is con- 
sidered to appear at this temperature. 

(iv) The dielectric constant increases below Tni = 37 
K. The lattice-Fcrri (T-?-|) state is considered to 

7 10 

appear.— 

(v) The neutron-scattering data of (1 1 1) begin to in- 
crease at Tn2 = 31 K. It is considered as the second 
magnetic phase transition as the temperature de- 
creases from the room temperature. The tempera- 
ture dependence of the (111) data is linear with 

(vi) All the neutron data saturate at T = 17 K, below 
which the magnetic order seems to be perfect*^ 

(vii) The dielectric constant also shows an anomaly at 
32 Kfi£ The temperature is very close to Tn2- It is 
not known whether it is another structural phase 
transition or not. 

The above experimental evidences are summarized in 
Fig.ffl 

The requirement (i) determines the energy scale of jf , 
the requirement (vi) determines jf , the requirement (iii) 
(Tni) determines the ratio jf /Jj 3 ) 12 ' 21 and the require- 
ment (ii) determines the ratio J\l J\- The other param- 
eters are determined by the temperature dependence of 
the structure factor between 20 K and 37 K. 



B. The spin-lattice model 

Numerical results are shown in Fig. [5] The seven pa- 
rameters are determined in order to fit the neutron data 
by visual inspection. The spin parameters are uniquely 
determined as 



-97 K, Jf 



-2.4 K, J 2 b = 0.14 K. 



(10) 



Those for the lattice system were not uniquely deter- 
mined. There are several choices that reproduce the ex- 
perimental results. We present two choices of the lat- 
tice parameters in this paper. Other possible parameter 
choices range between these two estimates. They are 



35 
30 
25 
20 
15 
10 
5 




x1(T 



(2/32mK 


Experiment 
Spin-Lattice — 












^1/3 i/3Tr~* 



































15 



20 



25 30 35 
Temperature [K] 



40 



CD 

CO 

'o 



15 
14 
13 
12 
11 
10 
9 
8 
7 
6 




Experiment ° 
Spin-Lattice — 

(b) 



J r L = 73 K, j\ 



-49 K. 



J, L 



0.38 K, A = 0.20, (11) 
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FIG. 5: Results of the spin-lattice model. The lattice param- 
eters are set to the values in Eq. (|11[) . (a) The MC data of the 
structure factor are compared with the neutron experimental 
data.^ The MC data of (ff) 2 are multiplied to coincide with 
the experimental data of (111). Those of (ff/ 3 ) 2 are multi- 
plied to coincide with the experimental data of (1/3 1/3 1) and 
(2/3 2/3 1). (b) The MC data of the uniform magnetic sus- 
ceptibility are compared with the experimental data.— The 
amplitudes of the simulation data and the constant contribu- 
tion from the non-magnetic impurity are determined so that 
the maximum value and the minimum value agree with the 
experimental data. 



and 



J f L = 61 K, J X L = -57 K, J 2 L = 0.61 K, A = 0.24. (12) 

Both parameter choices reproduce the neutron data and 
the susceptibility data fine. Quality of the fitting is the 
same within the visual inspection. These parameters re- 
produce the susceptibility data from 20 K to 140 K, in- 
cluding a convex change at 37 K. The intermediate phase 
between 31 K and 37 K is identified as the spin-PD phase, 
because the noncquilibrium relaxation data in Fig. [H] ex- 
hibits that ff disappears exponentially. 

We cannot uniquely determine the lattice parameters 
because of lack of information that determines . In 
the present lattice system, there is not an observable 
corresponding to the magnetic susceptibility, which de- 
termines jf . For each choice of j\, we can find j\ , J\, 
and A in order to satisfy the experimental data. The 
ratio Jc/Jc takes a value from 1.6 to 1.0, which agrees 
with the previous paper ^ 
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FIG. 6: A nonequilibrium relaxation plot of the magnetic 
structure factor /f when the simulation starts from the 
fcrrimagnetically-ordered state. The temperature is 35.5 K, 
just below Tni = 37 K. The relaxation function of the spin- 
only system converges to a finite value, while that of the spin- 
lattice system decays exponentially. The lattice parameters 
are those of Eq. (pH) . 



We have observed the sublattice magnetization and 
the sublattice polarization. Figures EJa) and [7]Jb) show 
the temperature dependence of the profiles. As shown 
in Fig. [TJb), the structural phase transition occurs at 
T ~ 90 K, below which the lattice system takes the 
lattice-PD (T-J.-0) state. It reproduces the experimental 
requirement (ii). The spin system remains paramagnetic 
in this temperature region. 

The spin transition and the lattice transition occur 
near 37 K (Fig. EJJa)). They are not always simultane- 
ous. When is set to 73 K, the spin transition occurs 
at 37.2 K, which is slightly higher than the lattice transi- 
tion temperature, 36.8 K. On the other hand, the lattice 
transition temperature becomes 38.5 K when is set 
to 61 K, while the spin transition temperature remains 
the same. The experimental finding of the simultaneous 
magneto-dielectric phase transition at 37 K is an acciden- 
tal coincidence. The two transition may occur at slightly 
different temperatures. 

In the case of J\ = 73K, a weak ferrimagnctic state 
appears at 37.2 K. This is because the lattice takes the 
lattice-PD (f-j-0) state at this temperature, and it favors 
the spin-Ferri state (Fig. [21(b)). When the lattice transi- 
tion occurs at 36.8 K, the spin-Ferri state disappears and 
the spin-PD state appears because the lattice-Ferri state 
favors the spin-PD state (Fig.f^a)). This is an outcome 
of the spin-lattice coupling. The situation changes when 
is set to 61 K. The lattice transition occurs at 38.5 
K. Since the lattice system takes the lattice-Ferri (T-t-J,) 
state at 37.2 K, the favored spin order is the spin-PD. 
The weak spin-Ferri state does not appear in this case. 
This weak ferrimagnctic phase is so narrow that it may 
not be observed in experiments if it exists. 

It is noticed that the spin transition temperature is ro- 
bust against the change of the lattice parameters. As far 
as we observed in the MC simulations with several lat- 
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FIG. 7: (Color online) Sublattice profiles of spin variables 
and lattice variables. The amplitude is normalized to unity 
when the sublattice order is perfect. Thin (red and green) 
lines are results of the = 73 K parameter set, while thick 
(magenta and light-blue) lines are results of the j\ — 61 K 
parameter set, respectively, (a) The temperature ranges from 
15 K to 40 K, where the corresponding neutron experiments 
were performed, (b) The high-temperature data. 



tice parameter choices, the spin transition temperature 
from the paramagnetic phase to the intermediate phase 
always occurs at 37 K. The spin transition temperature 
seems to be determined by the spin parameters Eq. (|10|) . 
On the other hand, the lattice system controls the type 
of the spin order at this temperature. 

In the intermediate phase, the amplitudes of two f- 
shiftcd sublattice polarization and that of one .[-shifted 
sublattice polarization are different. The former one is 
not saturated, while the latter is saturated. It is the 
two-sublattice lattice-Ferri (T-T--L) state (Fig. QJb)). An 
increase of the t-shifted polarization is slow and almost 
linear with the decreasing temperature. 

The low-tcmpcraturc magnetic transition occurs at 31 
K, below which the ferrimagnetic state appears. Each 
sublattice magnetization takes a different value. An in- 
version symmetry between a spin-up sublattice and a 
spin-down sublattice is broken. It is the thrcc-sublatticc 
ferrimagnetic state, which was predicted to appear by 
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FIG. 8: The structure factor data of the lattice system. The 
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the mean-field approximation^ As the temperature de- 
creases, the sublattice magnetizations approach the unity 
and the perfect ferrimagnetic order is realized at tem- 
peratures near T = 17 K. This saturation temperature 
depends on jf . 

The lattice transition always occur at the low- 
temperature magnetic transition temperature. It is the 
simultaneous spin-lattice phase transition. Above this 
spin-lattice transition temperature, the lattice system 
takes the two-sublattice lattice- Ferri state (Fig. Hlb)), 
which is the ground-state configuration. When the spin- 
Ferri state appears below the transition temperature, the 
lattice state is deformed toward the lattice-PD (|-|-0) 
state because the spin-Ferri state favors the lattice-PD 
state. A combination of the partial lattice- Ferri state and 
the partial lattice-PD state yields the three-sublatticc 
lattice-Ferri state (Fig. [ljd)). The space group changes 
from P3cl to P3cl. It is a clear evidence for a strong 
correlation between the spin system and the lattice sys- 
tem. The lattice system alone cannot make this struc- 
tural phase transition, because the lattice system is in 
the ground-state phase above the transition temperature. 
Therefore, this may be the spin-lattice transition driven 
by the spin degrees of freedom. 

Figure [8] shows the MC results of the structure factor 
of the lattice system. The /1/3-structure factor shows 
small anomalies at 37 K and 31 K. The temperature de- 
pendence between 37 K and 31 K is slow and almost 
linear with T . The /i-structure factor data qualitatively 
agree with the experimental results of the spontaneous 
polarization^ The data show a decrease below 31 K 
because the lattice-Ferri state is deformed toward the 
lattice-PD state by the ferrimagnetic transition. In the 
real experiment^ the spontaneous polarization shows a 
minimum at 23 K, while the minimum occurs at 28 K in 
the present simulation. 



C. The spin-only model 

We show that the spin-only system cannot explain all 
the experimental data. The three parameters, jf , jf , 
and jf , are determined in order to fit the neutron exper- 
imental data fine by visual inspection. Figure O shows 
the result. The estimates are 

J c s = -77 K, jf = -3.8 K, jf = 0.58 K. 

Agreement with the neutron data is good, while the 
susceptibility data disagree with the experiment signif- 
icantly. If we choose estimates that fit the susceptibility 
data, the neutron data disagree in turn. We cannot find 
estimates that satisfy both experimental data at the same 
time. 

The MC data of the structure factor of (1 1 1) take 
very small but finite values between 31 K and 37 K, as 
shown in Fig. [SJ It suggests that the ferrimagnetic order 
is finite and the intermediate spin-PD phase disappears. 
The direct transition from the paramagnetic phase to 
the ferrimagnetic phase is an outcome of a rather large 
estimate of Jf . 

The spin-only model does not consider the dielectric 
characteristics of RbCoBr3. We cannot explain the suc- 
cessive structural phase transitions by this model. Using 
three parameters, jf, jf, and jf, we only reproduce 
cither the neutron experimental data or the magnetic 
susceptibility data. On the other hand, the spin-lattice 
model explains quantitatively both dielectric and mag- 
netic properties using seven parameters. It suggests that 
the interplay between the spin system and the lattice 
system is essential in this compound. 



D. Perturbations 

We consider perturbation effects to the present model. 
It is intended to see how robust the characteristic be- 
haviors of RbCoBr3 are against perturbations as well as 
to propose further experimental investigations. Here, we 
consider three perturbations. The former two perturba- 
tions couple with the lattice system, while the last one 
changes the spin-lattice coupling parameter. 

First, a lattice interaction parameter is changed in or- 
der to make the lattice system hard. It corresponds to a 
pressure effect. We increase the interaction (spring con- 
stant) along the c-axis, J^, from 73 K to 97 K, while 
the other parameters remain the same as in Eq. pip . 
Figure [TU] shows the results. The lattice transition to 
the two-sublattice t _ T-| state occurs at 40 K (thin green 
lines in Fig. fTCTb)). while it occurs at 37.2 K (thin green 
line in Fig. [7{a)) in the unperturbed case. The magnetic 
transition occurs at 37.1 K (red and magenta lines in 
Fig. [TOT b)). and the PD state appears. This magnetic 
transition temperature is robust against the lattice per- 
turbation. The PD phase continues to the lower temper- 
ature, and the ferrimagnetic transition occurs at 24 K. 



9 



35 
30 
25 
20 
15 
10 
5 




xi(r 



(2/32/3l) ,,s S 


t 


Experiment ° 
Spin-only — 




















(1/3 1/3 ip- 





























15 20 25 30 35 
Temperature [K] 



40 



CD 

CO 

'o 



15 
14 
13 
12 
11 
10 
9 
8 
7 
6 




Experiment ° 
Spin-only — 

(b) 



20 40 60 80 100 120 140 
Temperature [K] 



FIG. 9: Results of the spin-only model, (a) The MC data 
of the structure factor are compared with the neutron exper- 
imental data.— The MC data of (/f ) 2 are multiplied to coin- 
cide with the experimental data of (1 1 1). Those of (/f/3) 2 
are multiplied to coincide with the experimental data of (1/3 
1/3 1) and (2/3 2/3 1). (b) The MC data of the uniform 
magnetic susceptibility are compared with the experimental 
data.— The amplitudes of the simulation data and the con- 
stant contribution from the non-magnetic impurity are deter- 
mined so that the maximum value and the minimum value 
agree with the experimental data. 



We observe the small bifurcation of the lattice profiles at 
this temperature. The two-sublattice lattice-Ferri state 
is deformed very weakly to the three-sublattice lattice- 
Ferri state. The simultaneous spin-lattice transition also 
occurs in the perturbed system. The PD phase becomes 
wider as in the typical ABX3 compounds. Another clear 
difference from the original parameter set is that the (1 
11) structure factor is convex when it appears at 24 K, 
while the original one is linear. 

Second, we changed J\ from 0.38 K to 0.75 K, while 
the other parameters remain the same as in Eq. pip . 
This perturbation favors the lattice-Ferri (T-T~l) state. 
It may correspond to applying the electric field to this 
compound. The results of the spin system are the same 
as the case where we changed J^. As shown in the fig- 
ures, both perturbations produce the same temperature 
dependences of the spin profiles. On the other hand, the 
lattice profiles are different. The lattice transition 
occurs at 45 K (thick light-blue lines in Fig.fTOTb')'). while 
it occurs at 40 K in the perturbation. Experiments 
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FIG. 10: (Color online) Perturbations on exchange interac- 
tions. In one case, we increased J\ from 73 K to 97 K, (thin, 
red and green) In the other case, we increased j\ from 0.38 
K to 0.75 K (thick, magenta and light- blue). The other pa- 
rameters are the same as in Eq. (|ll|l . (a) Neutron data of 
RbCoBr3 are compared with the MC data of the structure 
factor in the perturbed cases. The MC data are multiplied 
so that the saturation value coincide with the experimental 
data at 17 K. (b) Sublattice profiles of the spin and the lattice 
variables. 



under the electric field or the high pressure may detect 
these changes. 

The last perturbation changes the spin-lattice coupling 
parameter. We set A = 0.1 and 0.3 while the other pa- 
rameters are unchanged from Eq. (|11| . Figure fTTI shows 
the results of the structure factor compared with the neu- 
tron experimental data. The PD transition temperature 
and the saturation temperature are robust against this 
perturbation. On the other hand, the ferrimagnetic tran- 
sition temperature depends on A, which is observed by 
the change of slope of the (1 1 1) data. As we increase A, 
the ferrimagnetic transition temperature increases. The 
spin-lattice coupling relaxes frustration and stabilizes the 
ferrimagnetic state. 

Comparing the above with the results of the original 
parameters (Figs. [5] and [7]) we notice that there are sev- 
eral spin-lattice effects. The ferrimagnetic state appears 



10 



</) 35 

i, 30 

£ 25 

1 20 
o 15 
o 10 
5 5 

CD 

2 



xi(r 



(2/3 2/C 


ISs „^ 1 Experiment 
* Ir^T^^ A=0.3 - 
A=0.1 


(173 1/C 
















"(HI 







15 20 25 30 35 
Temperature [K] 



40 



FIG. 11: (Color online) Neutron data of RbCoBr3 are com- 
pared with the MC data of the structure factor when A is 
set to 0.1 and 0.3. The other parameters are the same as in 
Eq. (|lip . The MC data are multiplied so that the saturation 
value coincide with the experimental data at 17 K. 



at higher temperatures because of the relaxation of frus- 
tration. The (1 1 1) structure factor linearly depends 
on the temperature. There appears a three-sublatticc 
T-T-J. state. These characteristic behaviors in RbCoBr3 
are fragile and disappear when the spin-lattice coupling 
changes. Control of the lattice system by the electric 
field or the pressure may produce a new effect to the 
spin system. 



V. DISCUSSION 

The successive phase transitions of RbCoBr3 are well 
explained by the spin-lattice model introduced in this 
paper. Numerical data of our Monte Carlo simulations 
quantitatively agree with the experimental results. The 
spin-lattice coupling is found to be essential in this sys- 
tem. It produces nontrivial behaviors of RbCoBr3 differ- 
ent from other typical ABX3 compounds. The present 
analysis was enabled by the new cluster flip algorithm, 
which eliminates the slow MC dynamics in the quasi-one- 
dimensional frustrated spin system. 

The magneto-dielectric transition at 37 K is not always 
simultaneous. It is a coincidence that the spin transition 
and the lattice transition occur at the close temperatures 
in RbCoBr3. They may differ if the interaction parame- 
ters are different. The spin transition temperature is pos- 
sibly determined independently from the lattice system. 
The spin-lattice coupling only determines what type of 
the spin order is realized below this transition tempera- 
ture. 

On the other hand, the magneto-dielectric transition 



at 31 K is always simultaneous. It is the spin-driven lat- 
tice transition. The lattice symmetry changes in order to 
realize the ferrimagnetic state. Therefore, the anomaly of 
the spontaneous polarization observed experimentally 10 
at 32 K is considered as an indication of another struc- 
tural transition, where the space group changes from 
P3cl to P3cl. Further experiments to ensure this theo- 
retical prediction is expecting. An anomaly at 9 K ob- 
served experimentally has not been identified within the 
present spin-lattice model. 

The criticality of the phase transitions is an inter- 
esting future problem. The linear temperature depen- 
dence of the (1 1 1) structure factor below 31 K in 
Fig. OJa) may be an indication of the mean-field uni- 
versality j3 = 1/2. Though the interaction range is 
limited to the second-nearest neighbor, the spin-lattice 
coupling effectively makes it long-ranged because of the 
large correlation lengths of both spin variables and lat- 
tice variables. This mean-field criticality is supported by 
a model proposed recently by Miyashita et al^ Their 
model for a magnetic phase transition in spin-crossover 
materials is similar to our spin-lattice model. In their 
model, the spin takes either a high-spin state or a low- 
spin state. The volume of a magnetic ion depends on 
the spin state, which produces an effective spin-lattice 
coupling. They observed the mean-field universality by 
the detailed scaling analysis on the model system. If 
the mean- field universality appears in RbCoBr3, it may 
be observed in other magneto-dielectric compounds, e.g., 
i?Fe 2 4 ,2MI 

It should be commented that our mean-field-like treat- 
ment of the MC updating may have affected the critical 
phenomenon. This treatment averages the lattice vari- 
ables along the chain. It may produce effective long-range 
spin-spin correlations along the chain. 

In the present model, the lattice parameters have not 
been determined uniquely The lattice model is simple, 
taking only the elastic energy into account. Our assump- 
tion of the spin-lattice coupling only models the deforma- 
tion of the lattice system as an influence to the spin sys- 
tem. Some modifications to the model may be necessary 
when we discuss the magnetic-dielectric cross correlation 
under an electric field and a magnetic field. 
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